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ABSTRACT 


A  Blapl*  analysis  Is  presented  for  predicting 
pressure  lose  and  density  factors  for  two-phase, 
one-cosgionent  annular  flow  with  and  without  heat 
generation.  All  four  coahlnatlons  of  laminar  and 
turbulent  flow  In  both  the  annulus  and  core  are  con¬ 
sidered.  The  analysis  Is  based  on  assumed  velocity 
profiles  and  matching  velocity  and  shear  stress  at 
the  ll(|uld-vapor  Interface.  The  theory  Is  found  to 
comiiare  favorably  with  experimental  results  on  both 
vertldaX  waft  tiorlsontal  sections.  The  application 
of  the  thedry  to  natural  circulation  steam  genera¬ 
tors  Is  also  presented. 
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PREDICTION  OF  THE  PRESSURE  LOSS  AND  DENSITY 
FACTORS  FOR  TWO- PHASE  ANNULAR  FLOW  WITH 
OR  WITHOUT  KBAT  GENERATION^ 


By  J.  C.  Westmoreland 


NOMENCLATURE 

The  following  nomenclature  Is  used  In  the  paper: 

A  -  cross  sectional  area 
a  -  internal  radius 
D  -  Internal  dianeter 
F  -  sumnation  of  pressure  loss  factors 
f  -  moody  friction  factor  over  four 
g  -  gravity  constant 
h  -  enthalpy 

k  -  specific  pressure-loss  factor 
n  -  number  of  risers 
P  -  two-phase  pressure  gradient 
p  -  fluid  pressure 
q  -  heat  transfer  rate  per  imit  area 
r  -  volume  fraction  of  two-phase  constituent 
s  -  internal  flow  path  or  stream  line  coordinate 

V  -  velocity  of  fluid 

V  -  specific  volume 
w  -  weight  rate  of  flow 
X  -  vapor  qua3.ity 

d-  angle  of  inclination  from  vertical 

n  -  (1  -  V^) 

^  The  material  for  this  Paper  was  obtained  in  the  course  of  engineering  dewel- 
opnsnt  work  under  Contract  No.  W-31-109  Eng  52,  sponsored  by  the  General  Electric 
Company  which  operates  the  Knolls  Atomic  Power  Laboratory  for  the  Ikiited  States 
Atoudc  Energy  Comnission. 
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with  subscript  dynamic  viscosity,  otherwise  viscosity  ratio 
V  -  Vf  /Vg 
jO  -  fluid  density 
(S  -  water  level  swell  factor 

(1  -  Vf /Vg) 

'T -  shear  stress 

0  -  two-phase  flow  distribution  factor 
-  reclrctilatlon  ratio 

Subscripts 

d  -  dofwncomer  liquid 
f  -  Saturated  liquid 
g  -  saturated  vapor 

o  -  equivalent  liquid  flow  f ot  two-phase  flew 
w  -  solid  boundary  condition 


Note:  Average  values  of  velocity  are  denoted  by  a  cross  bar. 
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In  the  design  of  two-phase  flow  systems  for  a  wide  field  of  application, 
an  analytical  treatment  of  the  pressure  loss  and  density  variations  involved 
would  be  extremely  \i3eful;  in  some  cases  >here  experimental  results  are 
lacking  or  practical  experience  is  missing  this  treatment  becomes  quite 
necessary.  "Because  of  the  more  recent  en^jhasis  on  the  design  of  nuclear 
power  plants  which  comprise  liquid-cooled  reactors  and  steam  generators,  a 
greater  iiqportance  has  been  placed  on  the  need  for  such  a  treatment  in  the 
power  generation  field.  In  this  regard,  the  design  of  the  boiling  water 
reactors  and  nuclear  steam  generators  has  presented  problems  which  pertain  to 
their  control  and  performance.  The  solution  of  these  problems  requires  a 
detailed  knowledge  of  the  two-jhase  flow  phenomena  that  occur  in  the  associated 
steam  generation  sections.  Heretofore,  graphical  solutions  of  a  con^letely 
enq;>irical  nature  have  been  employed  exclusively  in  solving  these  particular 
problems;  however,  these  solutions  have  proven  iit^ractical  for  detailed 
analysis  to  be  made  with  an  extensive  variation  of  operating  conditions.  These 
techniques  offer  little  help  in  attempting  to  predict  the  transient  behavior 
of  any  specific  system. 


The  folloedng  analysis  of  two-phase  annular  flow  has  been  made  to  provide 
a  means  whereby  purely  analytical  procedures  may  be  en^loyed  In  predicting  the 
pressure  loss  and  density  relations  for  such  flows  in  systems  for  whldi  there 
Is  a  distinct  lack  of  experimental  data.  Throu^  the  use  of  establidied 
mechanics,  a  general  theory  for  two  phase  annular  flow  has  been  develop^  for 
evaluating  flow  problems  of  this  nature.  The  analysis,  however,  represents  a 
proposed  method  for  conrelating  two-phase  data  in  the  appropriate  range  of 
volime  fractions.  The  results  of  the  analysis  have  been  enployed  in 
successfully  predicting  the  experimental  data  of  a  number  investigations; 
the  more  general  results  have  been  reduced  to  dimensionless  form.  The 
practical  application  of  the  theory  to  a  specific  case  will  be  presented  and 
the  feasibility  of  the  results  is  being  studied  with  regard  to  the  design  of 
natural  circulation  steam  generators  for  nuclear  power  plants. 

TWO-mASE  FLCW  ANAUSIS 

The  basic  concept  behind  the  present  analysis  of  two-phase  annular  flow 
consists  of  prescribing  a  unique  criteria  in  terms  of  boundary  conditions  which 
give  the  recirculation  ratio  as  an  explicit  function  of  only  the  thermodynamic 
state  and  volume  fraction  in  lieu  of  the  particular  system  geometry  parameters 
and  absolute  flow  rate.  'Riis  important  result  is  obtained  by  performing  the 
following  integrations: 


where  Wf  is  the  liquid  annulus  mass  flow  rate  and  Wg  is  the  gas  core  mass  flew 
rate.  The  boundary  conditions  for  these  integrations  are  zero  liquid  velocity 
at  the  wall  and  equal  velocity  and  shoar  stress  at  the  two-phase  fluid  inter¬ 
face.  The  velocity  distribution  for  turbulent  flow  is  taken  to  be  the  <Mie- 
seventh  law,  and  for  laminar  flow,  the  square  law  is  enployed.  The  basic 
assunptions  for  this  analysis  may  be  enumerated  as  follows: 

!•  We  assume  constant  static  pressure  across  any  flow  cross  section. 

2.  We  assume  constant  fluid  properties  across  any  flow  cross  section. 

3.  We  neglect  the  effects  of  any  interphase  dispersion  or  entrainment. 

ii.  We  assume  saturated  generation  of  v^cr. 

It  has  been  recognized  that  each  phase  may  experience  a  transition  traa 
the  laminar  to  the  turbulent  flow  regime  and  vice  versa  under  different 
conditions,  and  the  indicated  integrations  have  been  made  for  all  permutations 
of  the  two  possible  simviltaneous  flow  conditions.  The  appropriate  velocity 
distribution  for  eadi  phase  has  been  selected  in  accordance  with  the  previously 
stated  criteria.  The  results^  of  these  integrations  are: 

Laminar  -  Laminar 


♦A  detailed  derivation  of  the  tirbulent-laminar  case  is  given  in  the  Appendix 


u 


Turbulent  -  Turbulent 


Laminar  •  Turbulent 

uhere 

Turbulent  -  Laminar 

where  the  two-phase  flow  distribution  factor  0  has  been  extracted  from  the 
result  obtained  for  mass  flow  ratio 
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In  any  one  permutation  of  flows  the  liquid-phase  flow  condition  is  described 
first.  The  flow  distribution  factor  is  dimensionless  and  independent  of  the 
system  geometry  parameters,  therefore,  it  becomes  only  a  general  function  of 
the  thermodjmamic  state  and  volume  fraction  of  gas.  VIhile  the  previous  derlva- 
tims  were  accomplished  for  a  tubular  cross  section,  the  results  are  felt  to 
be  representative  of,  and  have  been  centred  favorably  with,  the  flow  between 
two  flat  surfaces. 


The  flow  regime  associated  with  one  phase  or  the  other  may  be  ascertained 
once  we  are  able  to  predict  the  appropriate  Reynold's  number.  This  may  be 
done  for  the  gas  phase  taking  the  characteristic  velocity,  which  is 
descriptive  of  the  ratio  of  inertia  forces  to  the  viscous  forces  at  the  two- 
idiase  interface,  as  the  excess  of  average  gas  velocity  over  the  average  liqiiid 
velocity.  The  Reynold's  number  for  the  gas  phase  is  then: 

It  becomes  convenient  to  express  the  ratio  of  average  velocities  in  dimension¬ 
less  form;  this  ratio  may  be  obtained  from  the  results  of  Equation  (6)  as 
follows: 
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(8) 


For  the  liquid  phase  we  arbitrarily  define  the  Reynold's  number  in  terns  of  an 
apparent  liquid  flow.  These  terms  were  derived  by  integrating  the  velocity 
profile  in  the  annulus  across  the  entire  cross  section.  The  Reynold's  number 
for  the  liquid  phase  then  becomes 


(9) 


iDn  eocamlning  the  possibility  of  any  one  permutation  of  laainar  and 
txirbulent  flows  occurring  in  our  model,  a  useful  result  is  obtained  by  takii^ 
the  ratio  of  Equations  (9)  to  (7) 


For  any  one  thermodynamic  state  It  can  be  shown  that  at  moderate  voltw 
fractions  of  gas  (rg  >  0«U),  the  Reynold's  number  for  the  liquid  phase  is  much 
neater  than  that  for  the  gas  phase,  and  at  very  hi^  volume  fractions  of  gam 
(rg  ■  0«9)«  The  two-phase  flow  is  characterized  by  the  Reynold's  nuadier  for 
eadi  phase  becoming  approodmately  equal.  This  observation  would  indicate 
that,  without  prior  knowledge  of  the  geometry  car  flow  conditions  for  any  par¬ 
ticular  system,  the  most  probable  configuration  of  flow  regimes  would  be 
turbulent-laminar  at  moderate  volume  fractions  and  laminar-turbulent  at 
volume  fractions.  This  deduction  for  ^obable  mcxies  of  flow  has  been  borne 
out  by  the  results  obtained  by  Dengler^^'  on  a  particular  system.  In  these 
experiments,  photographs  of  tw -phase  annular  flow  were  taken  lAiich  clearly 
indicate  m  turbulent  annulus  and  laminar  core  at  volume  Aractlon  of  gas  of 
^;>proziniately  0,20$  the  observations  are  distinctly  reversed  for  the  photo¬ 
graph  ^own  at  higher  voliaoe  fractions  of  approximately  0,90, 

TWD-HASE  PRESSIEE  LOSS  ANiLTSlS 

The  pressure  loss  which  occurs  across  a  flow  section  with  two-phase  flow 
has  been  derived  by  superisqjosing  linearly  the  effects  caused  by  hrH  friction, 
net  mementum  flux.  In  a  vaq>Qr  generation  section  the  a8suBq)tian 
of  no  interphase  dispersion  is  subject  to  question,  and  it  is  obly  xeeognized 
that  (sonsideration  diould  be  given  to  the  change  In  momentum  that  arises  firom 
the  formation  of  vapor  bubbles  in  the  liquid  film  at  the  wall.  Throng  a 
vertical  section  the  static  pressure  is  assumed  ccwistant  across  a  flow  section, 
and  we  tacitly  e^ploy  the  same  assunq)tion  for  the  inclined  or  horizcmtal  ease. 
As  previously  stated,  the  fluid  properties  are  taken  as  constant  across  a 
cross  section. 

The  local  pressure  loss  caused  by  wall  friction  was  determinsd  by 
acknowledging  the  wall  shear  stress  ihlch  resulted  from  the  liquid  film  in 
terms  of  the  apparent  liquid  flow  as  foUowst 


(11) 
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In  t«nui  of  tho  total  liquid  ilow*  ve  prescribe  the  frieti<»al  pressure  loss 
or  wall  shear  stress  in  a  sinilar  nanner 


(12) 


the  total  liquid  flow  is  considered  at  the  sane  thermodynamic  state  as  the 
two-phase  flow*  dividing  the  two  previous  expressitms  we  obtain 


where  the  friction  factor  fgpp  is  evaluated  at  a  Reynold's  number  for  the 
iq}parent  liquid  flow  derived  previously*  We  may  proceed  to  write  the  general 
etxpresslon  for  the  Arlction  factor  in  terms  of  the  Reynold's  number  as 


If  C  is  an  empirical  constant  that  represents  the  friction  factor  variation 
for  the  liquid  flows  considered  idien  tte  appropriate  flow  regime  is  established, 
then  we  may  write  the  local  pressure-gradient  ratio  as 


By  inserting  in  this  expression  the  appropriate  Reynold's  number  for  the 
apparent  liquid  flow  along  with  the  appropriate  n\;^er  for  the  total  mass  flow 
in  the  liquid  state,  we  readily  derive  the  following  general  result  for  the 
turbulent  annulus 


(16) 


The  previous  expression  is  dimensionless  and  indo^wx^wn^  of  system  geometry 
and  total  flow  rates|  hence  by  selecting  o;>propriate  values  for  the  expment  m, 
we  may  proceed  to  obtain  dimensionless  plots  of  two-phase  friction  pressure 
loss  versus  the  volume  fraction  of  steam  at  constant  thermodynamic  states.  It 
has  been  found  ^at  an  eaqponent  of  l/U  represented  extremely  well  the  variation 
of  friction  factor  with  a  Reynold's  number  for  the  case  of  flow  alcxig  a  flat 
plate  in  the  turbulent  regime;  the  same  e3q>onent  has  been  found  satisfactory 
for  data  taken  from  turbulent  tube  flow*  We  have  in  a  laminar  flow  the  well- 
known  exponent  of  one*  It  is  interesting  to  note  that  if  we  assume  the 
Reynold's  noiber  for  the  annulus  flow  to  be  the  same  as  for  the  total  liquid 
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steam  flow 


Pig.  2 


Plg.l 


floWj  then  thie  exponent  becaaee  0,  axxi  all  possible  c^ses  may  be  described  with 
values  of  m  ranging  fron  0  to  !•  Without  knowledge  of  any  one  particxilar  system 
and  the  assodated  flow  conditions,  it  can  be  shown  that  an  exponent  of  l/k  gives 
a  fairly  accurate  estimate  of  the  two-phase  friction  pressure  loss*  This 
conclusion  w^s  established  from  the  results  shown  in  Figure  !• 

For  the  case  in  which  v^or  is  being  generated  alcuig  the  flow  channels,  a 
diange  in  volume  fraction  of  steam  is  experienced  from  inlet  to  outlet*  The 
functional  variation  of  volume  fraction  with  path  length  may  be  derived  from 
the  heat  flux  in  the  foUowlng  manner: 


The  volune  fraction,  rg,  at  any  distance,  s,  along  the  flow  path  is  then 
determined  from  the  total  flow,  Wq,  the  steam  flow  above,  and  the  correspondii^ 
value  of  Vg  In  the  appropriate  permutation  flow  regimes*  After  the 
distribution  of  rg  in  the  s  direction  has  been  obtained,  the  two-iiiase 
frictional  pre88urfr<drop  may  be  averaged  ever  the  path  length  by  the  following 
integration: 


(18) 
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This  result  is  then  eB|>loyed  to  predict  the  frictional  pressure  loss  for  heat 
generation* 

The  pressure  loss  attributed  to  a  change  in  monentum  flux  nay  be  derived 
for  an  incremental  secticm  by  uriting  the  following  expression  for  this  change 
idien  the  inlet  condition  is  zero  volume  fraction  of  vapor  and  the  outlet 
condition  is  ^ecifled  by  rgS 


n>ofrt. 


This  result  may  be  written  in  terms  of  a  multiplying  factor  R,  aS  follows: 


(19) 


where 


(20) 


(21) 


The  expression  R  is  dimensionless  and  with  the  previous  result  for  the  two- 
{^ase  mass  flow  ratio 


(22) 


it  becomes  an  e^qpllclt  function  of  volume  fraction  of  steam  and  thermodynamic 
state*  With  the  total  heat  generaticai  specified  then  the  dxsnge  in  momentum 
flinc  is  determined  from  the  exit  volume  faction  of  v£^or* 

The  pressure  gradient  required  to  balance  the  effect  of  gravity  may  be 
determined  for  the  general  case  of  heat  generation  and/or  variation  in  cross- 
sectional  area  by  siiqply  integrating  the  gravitational  component  betwera 
terminal  points 

- 37^ - ^ 

where  Cl  is  the  angle  of  inclination  from  the  vertical  and  is  a  known  calculable 
function  of  the  path-length  S »  This  result  is  governed  by  the  recirculation 
ratio  whidi  in  turn  fixes  the  volume  flraction  of  gas  at  any  section  for  the 
case  beldg  studied*  For  one-component  heat  generation  we  obtain  the  recircula¬ 
tion  ratio  at  any  cross  section  ftom 


9 


3^-  - 


(21*) 


troBi  lAiich  the  Tolume  fk'action  of  Taper  may  be  determined  throu^  the  results 
of  the  prerlous  two-^ase  flow  analysis. 

The  total  pressure  loss  for  two-phase  flow  is  then  obtained  by  sunning 
the  preylous  results  to  glye 

# 

c 

(25) 

c 


This  result  is  knowi  onco  the  ge<mietry  of  the  system  is  determined,  and  when 
appropriate,  the  heat  generation  characteristics  are  specified  aloi^  with  the 
total  liquid  flow.  To  predict  the  pressure  loss  associated  with  the  flow 
tturou^  yalves  and  fittings,  and  also  the  losses  in  the  entrance  and  exit 
regions  of  channels,  the  yelocity-head  conc^t  is  enyoked*  The  appropriate 
single-phase  pressure  loss  factor  being  applied  to  the  total  mass  flow  is 


SP. 


(/-yS/S)  J 


(26) 


where  K  is  the  single-phase  loss  factor, 

/>  - 

and 

_ 

0-/3^)/}  A 


(27) 

(28) 


In  toe  absence  of  a  better  concept  to  use  in  predicting  such  losses,  this 
result  seems  reasonable.  Perhaps  as  the  state  of  the  art  of  two-phaae  flow 
is  further  adyanced,  if>propriate  empirical  correlations  will  become  ayailable 
that  describe  this  pressure  loss  more  accurately. 

APPLICATICaS 

The  present  analysis  has  been  employed  effectiyely  in  analyzii^  the 
dynamics  of  natural  circulation  steam  generators  for  nuclear  power  plants.  A 
sdiematie  diagram  of  this  unit  is  shown  in  Figure  2  and  toe  major  con^onents 
are  riser,  downcomer,  eyaporator,  and  steam  drum.  The  circulation  of  fluid  is 
in  toe  direction  indicated  and  is  established  by  toe  generation  of  steam  in 
toe  eyaporator,  toich  in  turn  causes  the  grayity  term  in  the  riser  to  dlminleh 
and  so  creates  an  unbalance  with  toe  downcomer  grayity  term.  This  condition 
results  in  a  circulati(xi  of  magnitude  determined  by  the  frictlcsial  aixi  momentum 
losses  throu^out  the  flow  path.  The  water  leyel  maintained  in  the  steam  drum 
is  one  of  the  independexxt  yariables  along  with  the  magnitude  of  the  circulation. 
To  solye  this  problem  we  sum  the  pressure  losses  around  the  flow  circuit  and 
equate  the  result  to  zero 


TURBULENT  TURBULENT- 


W.g.3  Phase  flow  distribution  factor  H2O  at  500  psla 


HjO  AT  ^00  PSIA  ■ 


Pig. 5  Two-phase  slip  velocity  ratio  H2O 
at  500  pslE. 


Pig.  4 
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Ja  this  InTMtlgation  it  becones  neoessarj  o  obtain  knowledge  of  the  transient 
behavior  of  these  nnits*  It  was  foiund  that  the  solution  for  the  transient  case 
could  be  eaqployed  effeetirely’  for  steady-state  solutions  aljqply  by  placing  axyr 
reasonable  Yalue  of  the  independent  Tariable  (water  level)  into  the  transient 
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solution;  the  result  would  be  the  correct  value  for  the  steam  flow  and  thermo¬ 
dynamic  state.  The  latter  values  are  obtained  from  the  solution  to  the  energy 
equations  for  the  system. 

The  application  of  conservation  of  mementum  arounc'  the  steam  generator 


Pig. 6  Slip  velocity  ratio  prediction 
HgO  at  615  psla 

Exit  values  of  a  vertical  heat  generation 
section 


circuit  and  the  zero  stonmation  of  pressures  around  the  closed  loop  is  next 
considered.  The  result  becomes j  for  the  case  of  a  turbulent  riser  flow, 

(31) 

and  is  the  summation  of  two-iAiase  k  factors  for  the  riser  flow  circuit. 

In  taking  a  mass  balance  on  the  water  in  the  steam  drum  with  the  condition 
of  matched  feedwater  flow,  we  obtain 


(32) 
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The  two-phase  flow  in  the  riser  is  solved  by 


4^/  = 


The  solid  water  steam  drum  water  level  parameter,  is  obtained  fran  the 
volume  fraction  of  steam  in  the  risers  by 


(33) 


/.  -4  - 


(31.) 


Equations  (30),  (31),  (32),  and  (33)  have  been  combirjed  to  give  the  volvuue 
fraction  of  steam,  rg,  as  a  function  of  the  steam  pressiare  and  flew  rate,  and 
also  derivatives  thereof.  The  latter  variables  are  employed  as  known  parameters 
inasmuch  as  the  solutions  to  the  energy  equations  for  the  system  are  not  the 
subject  under  consideration.  With  these  terms  as  driving  functions,  the  finite 
difference  method  has  been  used  to  obtain  a  numerical  solution  of  the  resultant 
differential  equation  by  means  of  the  IBM-650  digital  computer. 


The  resxilts  of  the  previous  analysis  have  been  en^loyed  to  study  practical 
steam  generators  as  the  design  indicated;  the  conqparison  between  theory  and 
data  has  been  very  good. 


RESULTS 

The  two-phase  flow  distribution  factor,  0,  computed  for  the  anniilar  flow 
of  saturated  water  and  steam  at  500  psia,  is  shown  in  Figure  3  as  a  functi(»i 
of  the  volume  fraction  of  steam  rg.  The  four  possible  permutations  of  the 
laminar  and  the  twbulent  flow  regimes  are  Indicated,  and  it  is  interesting 
to  note  that  fer  the  laninar- turbulent  case  and  also  the  turbulent-laminar  case 
the  water  annulus  dominates  the  magnitude  of  0  at  low  volume  fractions  with  the 
steam  core  dominating  at  hi^  volume  fractions.  This  latter  observation  is 
made  when  the  same  flow  regime  persists  in  both  the  annulus  and  the  core.  As 
the  thermodynamic  state  influences  0  through  the  ratio  of  viscosities  for  steam 
and  water,  it  can  be  shown  that  this  effect  may  be  neglected  for  some 
applications  over  a  wide  range  of  pressure. 

The  two-phase  local  pressxire  loss  factor,  ,  computed  for  the 

annular  flow  of  saturated  water  and  steam  at  500  psia,  is  shown  in  Figure  U  as 
a  function  of  the  volume  fraction  of  steam,  r».  Again,  the  four  possible 
permutations  of  flow  regimes  are  indicated.  These  results  show  that  for  this 
situation  the  four  curves  are  quite  similar  up  to  a  volume  fraction  of 
approximately  0.70,  and  thereafter  the  laminar-laminar  case  and  the  turbulent- 
turbulent  case  become  divergent  in  opposite  directiois.  In  this  case  of  two- 
phase  fiction  losBfS^/S/S  >  th^nsodynamlc  state  has  a  strong  influence 
on  the  magnitude  of  the  results. 

The  two-{^a8e  slip  velocity  ratio  computed  for  the  previous  conditions  is 
shown  in  Figure  5  as  a  function  of  volume  fraction  of  steam.  The  discussion 
given  for  the  results  of  the  two-i^ase  flow  distribution  factor  0,  shown  in 
Figure  3>  applies  identically  to  this  factor  with  regard  to  its  behavior 
characteristics.  Also  shown  in  Figxu'e  6  is  the  e]i9)irical  curve  proposed  by 
Cook'2)  whldi  was  obtained  from  experimental  data  taken  in  a  vertical 
rectangular  section  (1/2  in.  wide,  2  in.  deep  and  2  ft.  long  -  approximate 
dimensions)  for  water  flowing  at  600  psi  with  uniform  heat  generatlcm.  The 
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▼olune  fraction  of  steam  was  detennined  by  the  radioactive  tracer  method  and 
the  correlation  is  reported  to  be  good  trithin  i  25^*  I*'  the  discussion  given 

previously  with  regard  to  the  relative  independence  of  slip  velocity  with 
density*  the  conqparisw  of  the  analytical  results  with  the  experimental 
results  at  600  psia  (the  viscosity  discrepancy  has  negligible  effects  in  this 
case)  is  considered  excellent  for  the  estimated  turbulent— laminar  and  turbulent- 
turb^ent  flow  regime  permutations. 

For  the  two-coB5)onent  flow  of  a  liquid  and  a  gas*  the  analytical  results 
have  been  ea^loyed  in  predicting  the  experimental  data.'*^'  In  these  experiments 
the  concurrent  flow  of  a  liquid-gas  mixture  in  a  pipe-line  contactor  was  studied 
flnri  both  oU-air  and  water-air  combinations  were  considered.  The  measurements 
taken  were  pressure-loss*  mass  flow  rates*  and  static  voltmie  fraction  of  gas. 
Photographs  of  the  two  phase  flow  were  taken  In  glass  flow  sections  in  order 
that  the  different  types  of  flow  could  be  observed. 

In  employing  the  analysis  to  predict  these  data  we  could  either  start  with 
ea^  of  the  two  mass  flow  rates  or  the  volume  fraction  being  known  and  then 
proceed  to  confute  the  remaining  variables.  A  conparison  of  the  analytical 
results  with  three  sets  of  experimental  data  is  presented  in  the  table  with 
the  volume  fraction  of  liquid  being  taken  as  the  known  variable  for  the  analysis. 
The  coaparison  is  considered  good  althou^  the  change  in  density  of  the  air  as 
it  passed  throu^^i  the  contactor  was  neglected  since  only  the  average  volume 
fractions  were  reported  for  the  experiments.  To  distinguish  between  the  four 
possible  permutations  of  flow  regies  the  Reynold's  numbers  based  on  the 
stperficial  velocities  were  used  as  a  rougix  go±dd*  It  was  reported  that 
annular  flow  occurred  for  the  experiments  at  volume  fractions  of  gas  from 
0.50  to  0.95  idxere  the  flow  path  occurred  principally  in  the  horizontal  plane. 

The  analytical  results  are  tentatively  restricted  to  volume  fracticais  of  gas 
tran  0.20  to  0.90  for  vertical  flow*  and  the  question  arising  from  the  assumpticxi 
of  constant  static  pressure  across  a  flow  section  woxild  leave  some  doubt  in 
regard  to  Uiis  range  being  applied  to  the  horizontal  case. 

In  general*  the  analytical  results  for  the  two-phase  pressure  loss  factor* 

*  have  been  enployed  in  conputing  this  term  versus  the  vapor  quality 
based  on  the  wel^t  rate  of  flow  at  pressures  of  li^.T*  500*  and  2500  psia. 

These  results  are  shown  in  the  curves  of  Figure  1  s\iperimposed  cm  the  same 
carve  proposed  by  Martinelll  and  Nelson. It  is  to  be  noted  that  the  analytical 
results  are  based  on  the  turbulent- laminar  pemiubation  of  flow  regimes.  A 
da^ed  line  is  drawn  across  these  curves  to  indicate  when  a  voltmie  fraction  of 
vapor  of  0.9  occurs  for  the  results  calcnilated  from  Equation  (B-20).  It  is 
tentatively  proposed  that  the  annular  flow  model  begins  to  deteriorate  in  this 
range.  The  agreenent  between  these  two  sets  of  curves  is  considered  mcjre  than 
fortuitous  ^  that  the  enpirical  curves  were  based  on  the  data  of  Martinelli 
and  Nelson the  condition  of  annular  flow  is  only  isplled  and  not  specifically 
stated.  It  is  not  surprising  to  note  that  at  pressures  around  one  atmosphere* 
the  analytical  results  are  applicable  only  at  very  low  vapcr  cpiality  based  on 
the  weij^t  rate  of  flow  althou^  the  static  volume  fraction  of  vapor  is  high 
at  about  0.9. 
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Prediction  of  Two-Congjonent  Two-Riase  Flow  Air-Water  Mixture 


Eiqperijttental  Values  _ Predicted  Valuee 
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0.593 

81.2 

the  apparatus  there  were  four  passes,  and  pressure  drop  measurements 
were  made  in  ea^  one.  Vhere  ilPp  is  the  two-phase  pressure  loss  and  the 
subscripts  Q  and  L  refer  to  the  gas  and  liquid  phase,  respectively,  and  X  is 
the  Martinelli  parameter. 


CONCLUSIONS 


The  success  of  the  annular  flow  theory  in  predicting  two-phase  flow 
I^otcmena  accurately  idiere  there  is  a  liquid  f  iln  aln^  the  wall,  but  the 
mid-stream  conditions  in  no  way  resemble  the  orderly  distribution  proposed  in 
the  previous  analysis,  nay  be  attributed  to  the  following  postulation.  For 
evory  case  of  two-phase  flow  where  the  confining  surface  is  wetted  by  the  liquid 
phase,  corresponding  hypothetical  two-phase  annular  flow  exists  whidi  is 
diaraeterised  by  a  uniform  annulus  of  liquid  and  a  core  of  gas  or  vs^or.  Further, 
the  appreciate  velocity  distribution  may  be  assigned  to  both  ahe  annulus  and 
the  core  with  the  boundary  conditions  at  the  two-phase  interface  b  eing  prescribed 
by  equal  velocity  and  shear  stress.  This  postulate  is  useful  when  we  apply  this 
analysis  to  the  case  of  two-phase  flow  in  the  entrance  region  of  a  noo-heat 
generation  section  where  we  would  not  normally  expect  a  fully  developed  velocity 
profile  to  exist  imtil  acme  distance  downstream.  Although  such  profiles  are 
utilised  in  the  analysis,  we  take  refuge  bdiind  the  previous  postulation  in 
that  the  derived  two-phase  distribution  factor  represents,  in  essence,  a  state¬ 
ment  of  the  condition  of  flow  continuity.  Althou^  the  flow  section  may  not 
be  long  enou^  to  afford  fully  developed  flow,  the  proper  volume  fraction  of 
gas  and  slip  velocity  ratio  will  be  given.  The  analytical  results  for  the 
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tvD-phase  friction  pressure  loss  are  subject  to  modification  in  this 
particular  case.  We  would  also  have  to  apply  similar  reasining  in  the  case 
of  horizontal  two-phase  flow  ^ere  because  of  the  effect  of  the  heavier 
liquid  annulus  the  film  thidcness  is  greater  at  the  bottom  of  the  section 
than  at  the  top.  Ih  summary^  the  two-phase  posulate  represents  a 
criterion  for  the  practical  liquid  film  two-iiiase  flow  in  that  it  may  be 
diaracterized  by  an  idealized  annular  flow  of  the  type  described. 

Figure  7  the  two-phase  flow  configurations  of  mixed  flow  anmQar 
flow  are  shown.  In  mixed  flow  two  distinct  regimes  may  be  characterized  as 
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Pig -7  Two-phase  flow  configurations 


follows:  for  volume  fractions  of  gas  less  than  0.20,  bubble  flew  occurs  witn 
bubbles  of  gas  unifc^rmly  dispersed  in  a  parent  llqiiid  flow;  for  volume 
fractions  of  gas  greater  than  0.90  fog  ^ow  occurs  with  siheroids  of  liquid 
imlformly  dispersed  in  a  parent  gas  flow.  Annular  flow  is  assumed  to  occur 
between  volume  fractions  of  gas  between  0.20  and  0.90,  and  althou^  it  is 
quite  possible  that  in  some  cases  a  permutation  of  bubble  flow  and  fog  flow 
will  occur  with  an  annulus  of  bubble  flow  and  a  CKare  of  fog  flow  the  effects 
of  interphase  dispersion  have  been  neglected  in  the  analysis.  As  the  volume 
fraction  of  gas  is  an  iBq)ortant  variable  in  the  study  of  two-phase  fleer,  it 
is  necessairy  to  distingui^  between  the  static  voltnne  fraction,  as  defimd 
here,  and  a  triMiiaT*  term  frequently  used  lAiidi  is  derived  from  the  ratio  of 
the  volmne  rates  of  flow.* 
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APPEND!!.  DHtIVATIOM  OF  TWO-fflASE  TURBULflU  LAMINAR  FLOW 

In  this  specific  case  the  liquid  annulos  is  taken  to  be  in  the  turbulent 
regime;  the  gas  core  is  laminar.  The  velocity  profile  for  the  liquid  phase 
is  taken  in  terms  of  the  y  coordinate  being  the  familiar  one-seventh  variation. 

Vf  •  (3-1) 


The  velocity  profile  for  the  gas  phase  is  taken  in  terms  of  the  r  coordinate 
in  accordance  with  the  square  law 

(»-2) 

The  total  liquid  mass  flow  is  determined  from 

(»-3)  ^ 

The  total  gas  flow  is  established  Srcax  a  similar  integration 


The  boundary  conditions  necessary  to  eliminate  the  maximum  v  elocity  terms  are 
established  by  taking  a^  In  terms  of  the  vdlume  fraction  of  gas 
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at  the  two-phase  interface 
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(B»5) 

(B^) 


and 


^  9y 


(a-7) 


The  coordinate,  a^,  i^ere  the  gas  velocity  goes  to  zero  is  also  elininated  by 
the  previous  boundary  conditions.  Ey  performing  the  indicated  integrations 
and  by  dividing  the  mass  flow  rates,  we  obtain 


£*6  .  J~ 

where 


(B-3) 


^  [  Hf-  7j  * 

The  total  flow  rate  is  then  siiqply  the  sum  of  the  gas  and  liquid  flow 

^  (B-io) 

Now,  for  the  local  pressxire  loss  caused  by  wall  fricticxi,  we  have  the  wall 
shear  stress  given 


(B-11) 


wliere  Vapp  ^  average  velocity  for  the  apparent  liquid  flow  derived  from 
integrating  the  velocity  profile  in  the  annulus  across  the  entire  cross  section. 
The  Arieticn  factor  (fapp)  is  based  on  the  Reynold's  nvnber  of  the  apparent 
liquid  flow  as  fdLlcwst 


'V 


(3-12) 


At  the  sane  section  we  define  a  local  pressiire  loss  which  represents  the 
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gradient  obtained  at  the  same  thermodynamic  state  with  an  equivalent  total 
liquid  flow  equal  in  magnitude  to  that  of  the  two-phase  flew.  In  the  turbulent 
regime  the  friction  factor  is  given  by 


(B-13) 


The  local  pressure  loss  at  any  section  is  then  given  by 


Sf^  -  (two-phase  pressure  loss)  (B-Ui) 

m  (equivalent  total  liquid  flow)  (B-15) 

^y  dividing  these  two  expressions  and  by  substituting  fer  the  appropriate 
frieticai  factor  we  obtain 


(B-16) 


By  substituting  for  the  appropriate  Reynold’s  number  we  obtain. 

s/3 


(B.17) 


The  apparent  liquid  velocity  is  obtained  from  the  following  integratiems : 


(equivalent  total  liquid  flow), 
-Jo 


(B-18) 


(B-19) 


By  performing  these  integrations  and  by  m^ing  the  appropriate  substitutions 
we  obtain 

( 


(B-20) 


T!’3  local  two-phase  pressure  loss  Is  then  given  by  this  eaqression  once  the 
reference  pressure  loss  is  determined. 


